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Abstract. According to the dynamical quasiparticle model (DQPM) - matched to reproduce 
lattice QCD results in thermodynamic limit, - the constituents of the strongly interacting quark- 
gluon plasma (sQGP) are massive and off-shell quasi-particles (quarks and gluons) with broad 
spectral functions. In order to address the electromagnetic radiation of the sQGP, we derive 
off-shell cross sections of — > 7* , qg — 5- 7* + g and qg — >■ j*q{qg 7*?) reactions taking into 
account the effective propagators for quarks and gluons from the DQPM. Dilepton production in 
In-|-In collisions at 158 AGeV is studied by implementing these processes into the parton-hadron- 
string dynamics (PHSD) transport approach. The microscopic PHSD transport approach 
describes the full evolution of the heavy-ion collision: from the dynamics of quasi-particles in 
the sQGP phase (when the local energy density is above ~ 1 GeV/fm'^) through hadronization 
and to the following hadron interactions and off-shell propagation after the hadronization. A 
comparison to the data of the NA60 Collaboration shows that the low mass dilepton spectra 
are well described by including a coUisional broadening of vector mesons, while the spectra in 
the intermediate mass range are dominated by off-shell quark-antiquark annihilation, quark 
Bremsstrahlung and gluon-Compton scattering in the nonperturbative QGP. In particular, 
the observed softening of the rriT spectra at intermediate masses (1 GeV < M < 3 GeV) is 
approximately reproduced. 



1. Introduction 

While ttie properties of hadrons are rather well known in free space (embedded in a 
nonperturbative QCD vacuum) the masses and lifetimes of hadrons in a baryonic and/or mesonic 
environment are subject of current research in order to achieve a better understanding of the 
strong interaction and the nature of confinement. A modification of vector mesons has been 
seen experimentally in the enhanced production of lepton pairs above known sources in nucleus- 
nucleus collisions at Super-Proton-Synchroton (SPS) energies [H [2]. This can be attributed 
to a shortening of the lifetime of the vector mesons p, uj and 4>. The question arises if the 
enhancement might (in part) be due to new radiative channels [3] from the strong Quark-Gluon 
Plasma (sQGP). The answer to this question is nontrivial due to the nonequilibrium nature of 
the heavy-ion reactions and covariant transport models have to be incorporated to disentangle 
the various sources that contribute to the final dilepton spectra seen experimentally. 



2. The PHSD approach 

To address the vector meson properties in a hot and dense medium - as created in heavy-ion 
collisions - we employ an up-to-date relativistic transport model, i.e. the Parton Hadron String 
Dynamics model |4J (PHSD) that incorporates the relevant off-shell dynamics of the vector 
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Figure 1. Diagrams contributing to the dilepton production from a QGP: (a) Drell- 
Yan mechanism, (b) gluon-Compton scattering (GCS), (c) vertex correction, (d) gluon 
Bremsstrahlung (NLODY), where virtual photons (wavy hues) spht into lepton pairs, spiral 
lines denote gluons, arrows denote quarks. In each diagram the time runs from left to right. 

mesons as well as the explicit partonic phase in the early hot and dense reaction region. PHSD 
consistently describes the full evolution of a relativistic heavy-ion collision from the initial hard 
scatterings and string formation through the dynamical deconfinement phase transition to the 
quark-gluon plasma (QGP) as well as hadronization and to the subsequent interactions in the 
hadronic phase. 

In the hadronic sector PHSD is equivalent to the Hadron-String-Dynamics (HSD) transport 
approach [5l El [7j that has been used for the description of pA and AA collisions from SIS to 
RHIC energies and has lead to a fair reproduction of hadron abundances, rapidity distributions 
and transverse momentum spectra. In particular, HSD incorporates off-shell dynamics for vector 
mesons - according to Refs. [8j - and a set of vector-meson spectral functions [10] that covers 
possible scenarios for their in-medium modification. The transport theoretical description of 
quarks and gluons in PHSD is based on a dynamical quasiparticle model for partons matched 
to reproduce lattice QCD results in thermodynamic equilibrium (DQPM). The transition from 
partonic to hadronic degrees of freedom is described by covariant transition rates for the fusion 
of quark-antiquark pairs to mesonic resonances or three quarks (antiquarks) to baryonic states. 

Various models predict that hadrons change in the (hot and dense) nuclear medium; in 
particular, a broadening of the spectral function or a mass shift of the vector mesons has been 
expected. Furthermore, QCD sum rules indicated that a mass shift may lead to a broadening 
and vice versa [9J; therefore both modifications should be studied simultaneously. Thus we 
explore three possible scenarios: (1) a broadening of the p spectral function, (2) a mass shift, 
and (3) a broadening plus a mass shift. The HSD (PHSD) off-shell transport approach allows to 
investigate in a consistent way the different scenarios for the modification of vector mesons in a 
hot and dense medium. In off-shell transport the hadron spectral functions change dynamically 
during the propagation through the medium and evolve towards the on-shell spectral function 
in the vacuum. 

As demonstrated in Ref. [TO] the off-shell dynamics is particularly important for resonances 
with a rather long lifetime in vacuum but strongly decreasing lifetime in the nuclear medium 
(especially uj and </> mesons) but also proves vital for the correct description of dilepton decays 
of p mesons with masses close to the two pion decay threshold. For a detailed description of the 
various hadronic channels included for dilepton production as well as the off-shell dynamics we 
refer the reader to Refs. [TOllII]. 

3. SPS energies 

By employing the HSD approach to the low mass dilepton production in relativistic heavy-ion 
collisions, it was shown in that the NA60 Collaboration data for the invariant mass spectra 
for p'^p^ pairs from In+In collisions at 158 A-GeV favored the 'melting p' scenario [13J. Also 




Figure 2. The HSD results for the mass differential 
dilepton spectra from In + In collisions at 158 
A-GeV in comparison to the excess mass spectrum 
from NA60 [l3]. The actual NA60 acceptance filter 
and mass resolution have been incorporated [23]. 
The solid lines show the HSD results for a scenario 
including the collisional bradening of the p-meson 
whereas the dashed lines correspond to calculations 
with a 'free' p spectral functions for reference. The 
dash-dotted lines represent the HSD calculations 
for the 'dropping mass + collisional broadening' 
model. The bands represent the preliminary PHSD 
results incorporating direct dilepton radiation from 
the QGP in addition to a broadened /9-meson. 



Figure 3. Acceptance corrected mass 
spectra of the excess dimuons from In + 
In at 158 AGeV integrated over px in 
0.2 < pt < 2.4 GeV from PHSD compared 
to the data of NA60 [15j. The green 
dash-dotted line shows the dilepton yield 
from the vacuum p meson. The blue 
dashed line is the contribution to the 
dilepton yield from the in-medium p with 
broadened spectral function. The red solid 
line presents the sum of the in-medium p 
and QGP dilepton radiation (the latter is 
calculated in the on-shell approximation in 
this work). 



the data from the CERES Collaboration p[l] showed a preference for the 'melting p' picture. 
On the other hand, the dilepton spectrum from In+In collisions at 158 A-GeV for M > 1 GeV 
could not be accounted for by the known hadronic sources (see Fig. 2 of pT]). Also, hadronic 
models do not reproduce the softening of the niT distribution of dileptons for M > 1 GeV [13j. 
This observation pointed towards a partonic origin. 

Dilepton radiation by the constituents of the strongly interacting QGP proceeds via the 
elementary processes illustrated in Fig. [1} the basic Drell-Yan q + q annihilation mechanism, 
Gluon Compton scattering {q + g — )■ 7* + (7 and q + g — )■ J* + q), and quark -|- anti-quark 
annihilation with gluon Bremsstrahlung in the final state (q + q ^ g + J*)- In the on-shell 
approximation, one uses perturbative QCD cross sections for the processes listed above. The 
obtained off-shell elementary cross sections have been implemented into the PHSD transport 
code. Note that in our calculations the running coupling as depends on the local energy density 
e according to the DPQM [T2], while the energy density is related to the temperature T by the 
IQCD equation of state. Numerically, we observe that as is of the order 0(1) and thus the 
contribution of the higher order diagrams in Fig. [1] is not subleading! 

As we find in Fig. [3] the NA60 data favor the scenario of the in- medium broadening of vector 
mesons. Note that in the data - presented in this plot - the D-meson contribution has not been 
subtracted. The NA60 collaboration has published acceptance corrected data with subtracted 
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Figure 4. Left: Transverse mass spectra of dileptons for In+In at 158 A-GeV. Right: Right: 
The inverse slope parameter T^ff of the dimuon yield from In+In at 158 A-GeV as a function of 
the dimuon invariant mass in PHSD compared to the data of the NA60 Collaboration [1'6\ I15j. 

charm contribution recently [T5]. In Fig. [3] we present PHSD results for the dilepton spectrum 
excess over the known hadronic sources as produced in In + In reactions at 158 AGeV compared 
to the acceptance corrected data. We find that the spectrum at invariant masses below 1 GeV 
is well reproduced by the p meson yield, if a broadening of the meson spectral function in the 
medium is assumed. On the other hand, the spectrum at M > 1 GeV is shown to be dominated 
by the partonic sources. 

It is also interesting to note that accounting for partonic dilepton sources allows to reproduce 
in PHSD (cf. Fig. [3l rhs) the intriguing finding of the NA60 Collaboration [13l [15] that the 
effective temperature of the dileptons (slope parameters) in the intermediate mass range is 
lower than the Tefj of the dileptons from the hadronic phase. The softening of the transverse 
mass spectrum with growing invariant mass reflects that the partonic channels occur dominantly 
before the collective radial flow has developed. This feature of the data is also reproduced in 
PHSD (cf. Fig. m rhs). A detailed look at the PHSD results shows that in total we stih 
underestimate the slope parameter T^ff which might be due to missing partonic initial state 
effects or an underestimation of partonic flow in the initial phase of the reaction. 

4. RHIC energies 

In 2008, the PHENIX Collaboration has presented first dilepton data from pp and Au + Au 
collisions at Relativistic-Heavy-Ion-Collider (RHIC) energies of ^/s = 200 GeV [I?] which show 
an even larger enhancement in Au + Au reactions (relative to pp collisions) in the invariant mass 
regime from 0.15 to 0.6 GeV than the data at SPS energies |13[ll4j. We recall that HSD provides 
a reasonable description of hadron production in Au + Au collisions at ^/s = 200 GeV [19] such 
that we can directly continue with the results for e~^e~ pairs which are shown in Fig. [5] for p + p 
collisions. We find that the dilepton radiation in the elementary channel is well under control 
in HSD, which is equivalent to PHSD for p + p reactions. 

We step on with the case of inclusive Au + Au collisions in comparison to the data from 
PHENIX [T7J (Fig [6]). When including the in- medium modification scenarios for the vector 
mesons, we achieve a sum spectrum which is only slightly enhanced compared to the 'free' 
scenario. Whereas the total yield is quite well described in the region of the pion Dalitz decay 
as well as around the u and (j) mass, HSD clearly underestimates the measured spectra in the 
regime from 0.2 to 0.6 GeV by an average factor of ~ 3. The low mass dilepton spectra from 
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Figure 5. The HSD results for the mass 
differential dilepton spectra in case of pp 
collisions at ^/s = 200 GeV in comparison 
to the data from PHENIX [T^. The actual 
PHENIX acceptance and mass resolution have 
been incorporated (see legend for the different 
color coding of the individual channels) . Figure 
taken from 
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Figure 6. The HSD results for the mass differential dilepton spectra in case of inclusive 
Au + Au collisions at ^/s = 200 GeV in comparison to the data from PHENIX [T7]. The 
actual PHENIX acceptance filter and mass resolution have been incorporated [18J. In the upper 
part the results are shown for vacuum spectral functions (for p,u},<j3) including the channel 
decompositions (see legend for the different color coding of the individual channels). The lower 
part shows a comparison for the total e~^e~ mass spectrum in case of the 'free' scenario (dashed 
line), the 'collisional broadening' picture (solid line) as well as the 'dropping mass + collisional 
broadening' model (dash-dotted line). Figure taken from 



Au + Au collisions at RHIC (from the PHENIX Collaboration) are clearly underestimated in the 
invariant mass range from 0.2 to 0.6 GeV in the 'collisional broadening' scenario as well as in the 
'dropping mass + collisional broadening' model, i.e. when assuming a shift of the vector meson 
mass poles with the baryon density. We mention that our results for the low mass dileptons 
are very close to the calculated spectra from van Hees and Rapp |20j as well as Dusling and 
Zahed [21] (cf. the comparison in Ref. [22j). 

Consequently we attribute this additional low mass enhancement seen by PHENIX to non- 
hadronic sources. This assumption has to be tested quantitatively. However, in order to make 
quantitative predictions at experimentally relevant low dilepton mass (< 0.6 GeV) and strong 
coupling (05 ~ 0.5-^1), we have to take into account not only the higher order pQCD reaction 
mechanisms, but also the non-perturbative spectral functions and self-energies of the quarks. 
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Figure 7. Dilepton production cross sections 
in the Drell-Yan channel {q + q ^ fi'^ + fi~). 
The short dashes (black) line shows the on- 
shell, i.e. the standard pQCD result. The 
other lines show the off-shell cross section, in 
which the annihilating quark and antiquark 
have finite masses mi and rn-2 with different 
values: mi = 0.3 GeV, m2 = 0.15 GeV (solid 
magenta line), mi = 0.3 GeV, m2 = 0.3 GeV 
(dash-dotted red line), mi = 0.3 GeV, 772-2 = 
0.6 GeV (dashed blue line), mi = 0.6 GeV, 
m2 = 0.6 GeV (dash-dot-dot green line). 



Figure 8. Comparison of off-shell and on- 
shell cross sections for dilepton production in 
the gluon Bremsstrahlung q + q ^ g + 
channel at ^/s = 4 GeV. The dashed black 
line shows the on-shell (pQCD) cross section 
regularized by the cut on the gluon mass ficut = 
0.206 GeV, the blue dash-dotted line presents 
the off-shell cross section for the gluon mass 
fixed to /i = 0.8 GeV and on-shell quark and 
anti-quark (mi = 77x2 = m^ = 0). The red solid 
line gives the off-shell result for /i = 0.8 GeV, 
mi = 7712 = 777,3 = mg = 0.6 GeV. 



anti-quarks and gluons thus going beyond the leading twist. 

For this purpose, off-shell cross sections were derived in [24j for dilepton production in the 
reactions q + q ^ (Drell-Yan mechanism), q + q ^ g + (quark annihilation with the 
gluon Bremsstrahlung in the final state), q{q) + 5 — )• q{q) + l^l^ (gluon Compton scattering), 
g — )• q+q+l^l^ and q{q) — )• q{q)+g+l^l~ (virtual gluon decay, virtual quark decay) in the sQGP 
in effective perturbation theory by dressing the quark and gluon lines in the diagrams in Fig. 1 
with the DQPM propagators for quarks and gluons. The DQPM describes QCD properties in 
terms of single-particle Green's functions (in the sense of a two-particle irreducible approach) 
and leads to the notion of the constituents of the sQGP being effective quasiparticles, which are 
massive and have broad spectral functions (due to large interaction rates). 

Dilepton production cross sections in the Drell-Yan mechanism are plotted on the l.h.s. 
of Fig. El The short dashes (black) line shows the on-shell, i.e. the standard perturbative 
result. The other lines show the off-shell cross section, in which the annihilating quark and 
antiquark have finite masses mi and 7712 with different values: mi = 0.3 GeV, r?72 = 0.15 GeV 
(solid magenta line), 777i = 0.3 GeV, 7772 = 0.3 GeV (dash-dotted red line), 7771 = 0.3 GeV, 
7712 = 0.6 GeV (dashed blue line), r?7i = 0.6 GeV, 7772 = 0.6 GeV (dash-dot-dot green line). 

In Fig. [8l the off-shell cross sections for the quark annihilation with gluon bremsstrahlung 
channel at various values of quark and gluon off-shellnesses (masses) are compared to the on-shell 
(pQCD) result. The dashed black line shows the on-shell cross section regularized by a cut on 
gluon mass ficut = 0.206 GeV, the red solid line presents the off-shell cross section for the gluon 
mass fixed to // = 0.8 GeV and on-shell quark and anti-quark (mi = 7772 = 777,3 = 0). The blue 
dash-dotted line gives the off-shell result for = 0.8 GeV, 777,1 = "7,2 = m^ = mq = 0.6 GeV. 



One readily notices the shift of the maximum allowed mass of the lepton pair to a lower value 
(in order to produce a massive gluon in the final state). 

By implementing the off-shell partonic processes into the PHSD transport approach, a 
consistent calculation of the dilepton production in heavy-ion collisions at RHIC energies will 
be performed in near future. The comparison to the dilepton data double differentially in mass 
and pt will open the possibility to study the relative importance of different processes in the 
dilepton production and guide us towards a better understanding of the properties of matter at 
high densities and temperatures as created in heavy-ion collisions. 

5. Summary 

Dilepton production in relativistic heavy-ion collisions is a valuable probe of two phenomena 
of fundamental interest: the vector meson properties in a hot and dense medium and of the 
quark and gluon dynamics in the deonfined state of matter (sQGP). The Parton Hadron String 
Dynamics (PHSD) transport approach incorporates the relevant off-shell dynamics of the 
vector mesons as well as the explicit partonic phase in the early hot and dense reaction region. 

A comparison of the transport calculations to the data of the CERES and NA60 
Collaborations points towards 'melting' modification of the p'-meson at high densities, i.e. 
the broadening of the vector meson's spectral function. On the other hand, the spectrum for 
M > 1 GeV is shown to be dominated by the partonic sources. 

The low mass dilepton spectra from Au + Au collisions at RHIC (from the PHENIX 
Collaboration) are clearly underestimated by the hadronic channels in the invariant mass range 
from 0.2 to 0.6 GeV in the 'collisional broadening' scenario as well as in the 'dropping mass + 
collisional broadening' model, i.e. when assuming a shift of the vector meson mass poles with 
the baryon density. We attribute this additional low mass enhancement seen by PHENIX to 
radiation from the strongly interacting QGP. This assumption has to be tested quantitatively. 
Strong interaction of partons - reflected in their high masses and broad widths - lead to higher- 
twist corrections to the perturbative cross sections especially at the experimentally relevant low 
dilepton masses (< 0.6 GeV). By implementing the off-shell partonic processes into the PHSD 
transport approach, a consistent calculation of the dilepton production in heavy-ion collisions at 
RHIC energies will be performed in near future that will allow us answer whether the observed 
access can be attributed to the radiation by off-shell partons. 
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